Polarization probes of vorticity in heavy ion collisions 
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We discuss the information that can be deduced from a measurement of hadron (hyperon or vector 
meson) polarization in ultrarelativistic nuclear collisions. We describe the sensitivity of polarization 
to initial conditions, hydrodynamic evolution and mean free path, and find that the polarization 
observable is sensitive to all details and stages of the system's evolution. We suggest that an 
experimental investigation covering production plane and reaction plane polarizations, as well as 
the polarization of jet-associated particles in the plane defined by the jet and particle direction, can 
help in disentangling the factors contributing to this observable. Scans of polarization in energy and 
rapidity might also point to a change in the system's properties. 

PACS numbers: 25.75.-q, 13.88. +e, 12.38.Mh, 25.75.Nq 
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Parity violation, together with the self-analyzing na- 
ture of hyperon decays, provides us with the opportunity 
to study the polarization of hyperons produced in heavy 
ion collisions: In the rest frame of the hyperon Y, the 
angular decay distribution w.r.t. the polarization plane 

is [J 



dN 



1 + UyPy COS(0) 



(1) 



where ay is a hyperon-specific constant (measured in 
elementary processes [l[ ) , Py is the hyperon polarization 
and is the angle between the proton momentum and A 
polar axis. Similar relations arise, without the need for 
parity violation, in the strong decays of vector mesons to 
"p-wave" final states Q. 

If a specific net (Py) ^ exists in any axis definable 
event-by-event, it is in principle possible to measure it 
using Eq. Q] and the observed spectra of A, EE, f2 decay 
products. This opens a new avenue to investigate heavy 
ion collisions, which has been proposed both as a signal 
of a deconfined regime [1,0,-5, 6] and as a mark of global 
properties of the event [7|, [3, la, [lCj . 

Since QCD contains a spin-orbit coupling, a non-zero 
hyperon polarization in direction i is in principle present 
whenever the angular momentum density in that direc- 
tion, (l x T 0t ^ , is non-zero (T^ is here the energy mo- 
mentum tensor). As we will see, in both elementary p-p, 
p-A and A-A collisions is is possible to define directions 
where this vector might have non-zero components. 

The potential of hyperon polarization as a signal for 
deconfincment comes from the strong transverse polar- 
ization of hyperons in the production plane (left panel of 
fig. [ljobserved in (unpolarized)p + p and p + A colli- 
sions|llT [l2| . As suggested in d,|4, 5i], the disappearance 
of this polarization (which we shall call Py) could sig- 
nal the onset of an isotropized system where, locally, no 
reference frame is preferred. Something close to what is 



today called a "Quark Gluon liquid" (strongly interacting 
Quark-Gluon Plasma, sQGP). So far, no such measure- 
ment exists at RIIIC energies, through an Alternating 
Gradient Synchrotron (AGS) measurement [l3| yields a 
negative result (transverse polarization is comparable to 
p-A collisions). 

It has also been suggested 0] to use hyperon polar- 
ization in the reaction plane (left panel, Fig. [I} to test 
for local angular momentum in the matter produced in 
heavy ion collisions. The idea is that the initial momen- 
tum gradient in non-central collisions should result in 
a net angular momentum (shear) in this direction, that 
will be transferred to hyperon spin via spin-orbit cou- 
pling (This polarization direction will be called Py). A 
similar, through quantitatively different result, can be 
obtained from a microcanonical ensemble with a net an- 
gular momentum 

The STAR collaboration has recently measured the re- 
action polarization [14], reporting results consistent with 
zero. The production plane polarization measurement is 
also planned. 

In this paper, we make a few general considerations 
regarding the insights that can be gained from polariza- 
tion measurements. We examine how the polarization, in 
both production and reaction plane, is sensitive to initial 
conditions, hydrodynamic evolution, and mean free path. 
We suggest measuring polarization in different directions 
(production, reaction, and jet axis in jetty events) could 
provide a way to go beyond model-dependence. 

While throughout this paper we use the A polarization 
as our signature of choice, the points made here can be 
easily generalized to the detection of polarization of vec- 
tor mesons, also used as probes of polarization in a very 
similar way to hyperons [8( 
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FIG. 1: (color online) Definition of production and reaction plane. The beam line (traditionally the z axis) is perpendicular to 
the sheet. The dotted line, with arrow, indicates the direction of polarization of the produced A. 



INITIAL CONDITIONS AND REACTION PLANE 
POLARIZATION 

The QCD spin orbit coupling is capable of transform- 
ing the orbital angular momentum density [x x T®' 1 ^ — 

(x x p) into spin. The fact that the relevant quantity 
is angular momentum density, rather than its absolute 
value, can be seen intuitively by the requirement of lo- 
cality. Formally, it is apparent if the polarization from 
scattering is calculated explicitly using a Wigner Func- 
tion formalism 

For a large system, such as a heavy nucleus, we have to 
convolute the net polarizing interaction cross-section per 
unit of transverse nuclear surface (dAa/d 2 x±, where x± 
are the two directions perpendicular to the beam axis) 
calculated in [7| with the (initial) parton phase space 
distributions f(x±,p) to obtain the net local polarized 
parton phase space density p P R produced in the first in- 
teractions 

P P r(x±,p) = J d 2 x' ± d s p'f(x± - x' ± ,j) J^r 3 ') 

(2) 

where f(x±,p) is the local parton distribution of the 
medium. 

Provided the initial Debye mass and constituent quark 
mass is small, the quark polarization in the reaction plane 

(If) = J d 2 x ± d 3 pp PqR (x ± ,p) 

becomes Q 



p(x±) = / d 3 pf(x±,p) 



where 



is the participant transverse density and ft is a unit vec- 
tor perpendicular to both x± and p. In ultrarelativistic 
collisions all significant initial momentum is in the beam 
(z) direction. 

In non-central collisions with a finite impact parameter 
b, (p x x±) oc b ^ 0, thereby generating a resulting net 
polarization. 

Thus, the initially generated amount of reaction plane 
polarization is strongly dependent on the initial density- 
momentum correlation within the system. In other 
words, the reaction plane polarization could be a use- 
ful signature for probing the initial conditions within the 
system created in heavy ion collisions. 

According to the Glauber model, the initial density 
transverse coordinate distribution is given by the sum of 
the participant and target density pp,r- 



where 



Pp,t = T p ,t x± =F 



p(x±) = (pp(xj_) + Pt(x±)) <l>{y, rj) 



1 — exp 



(4) 



-ctnT t , p x± ± 



(5) 

and (Tat, Tp.r and b refer, respectively, to the nucleon- 
nucleon cross-section, the nuclear (projectile and target) 
density, and the impact parameter. 

How this density is longitudinally distributed in space- 
time rapidity 



(P q R ) 



d 2 x±p(x±)p.(xl X ft) ~ - (p z xj_) (3) 
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and flow rapidity 
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(the form of (f)(y, T])) is a crucially important model pa- 
rameter. The calculation of the hyperon polarization in 
the reaction plane 0] is dependent on an assumption of 
initial condition described in Fig. 1 and Eq. 1 of [7(. 
Such an initial condition (generally referred to as the 
"firestreak model"), is roughly equivalent to two "pan- 
cakes" , inhomogeneous in the transverse coordinate x± , 
locally inelastically sticking together. Each element of 
this system then streams in the direction of the local net 
momentum (Fig. [21 right column). 

Since projectile and target have opposite momenta in 
the center of mass frame, assuming projectile and target 
nuclei to be identical yields 



(j>(y, rj) ~ S(r))S (y - y cm (x±)) 



(8) 



where y cm is the local (in transverse space) longitudinal 
rapidity, corresponding to the flow velocity v cm given by 
momentum conservation. Thus 



<P,*x> ~ ^ {D P ) 
mjv 



where 



(D p ) = J d 2 xj_xj_ [pp(x±) 



(9) 



- Pt{x ± )] (10) 



If the colliding nucleons are also the interacting degrees 
of freedom, then the constant would be independent of 
energy. However, we know that at high energy the phys- 
ical degrees of freedom are partons, and the amount of 
partons each nucleus "fragments" into is highly energy 
dependent. 0] takes this into account using the energy- 
dependent parameter c(s). If one assumes all entropy to 
be created in the initial moment, c(s) as a function of 
energy can be estimated from final multiplicity using the 
well-known phenomenological formulae [15j | (~ 18 at top 
RHIC energies). 
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Assuming that all partons receive an equal share of mo- 
mentum, we get 



(PzX±) 



c(s)mjv 



(12) 



Since all nuclei have the same y/s, (PF) should be finite 
and constant over rapidity (Fig. 1 of pj). 

The physical validity of such a picture is compelling 
at low energies, when the baryon stopping of nuclear 
matter is large. At high energies and initial transparen- 



cies, however, a more generally accepted ansatz for initial 
condition is that approximated by a Brodsky-Gunion- 
Kuhn (BGK) [l(| [lTj picture, where the initial partons 
are produced ail-throughout the longitudinal flow rapid- 
ity spanned between the forward-traveling projectile and 
the backward-traveling target (middle panel of Fig. [2]). 
The spacetime rapidity 

is, in this picture, equal to the flow rapidity (Hub- 
blc/Bjorken expansion). 

If pp(x±) — Pt(x±), this reduces to a boost-invariant 
initial condition. For a non-central collision, however, 
such equality will only hold at the midpoint in x± of the 
collision region. Interpolating linearly in the rapidity y 
between pp (at y — y^) and p^ (at y = — j/i), we have 



4>{y,r}) = (A + yB)6(y-r)) 



Pp{x±) - Pt(xj_) 1 



pp(x±) + p T x± 2y L ' 
In particular, it means that for identical nuclei 



(13) 



(14) 



{PzX. 



dy smh(y)x±p(x±,y)dxL 
« (Dp) (y L cosh(y L ) - sinh(y L )) 



(15) 



(Note that no c(s) is necessary here since we are dealing 
with collective longitudinal flow). For this initial condi- 
tion, the axial symmetry of the initial pancakes forces 
the net polarization at mid-rapidity to be zero, and to 
rise as 0(yf j ) for y^ < 1. 

The rapidity distributions are summarized in Fig. [21 
and the corresponding shear created is summarized in (a) 
and (b) of Fig. H 

Thus, at high (RHIC and LHC) energies we expect net 
polarization around the reaction plane of A-A collisions 
should vanish at mid-rapidity and reappear in the target 
and projectile regions (Panel (a) of Fig. [3]). At lower 
energies, on the other hand, reaction plane A or A po- 
larization should be more uniform in rapidity space, and 
be significantly above zero at mid-rapidity (Panel (b) of 
Fig. EJ. 

Realistic nuclear geometries should not alter these very 
basic considerations, although for the BGK case they 
might considerably slow down the shear rise in rapidity. 
This is also true for corrections to linear interpolation in 
rapidity space. Detailed hydrodynamic simulations (l8j 
also reinforce the conclusion that within the boost invari- 
ant limit vorticity is negligible. 

A net vorticity will reappear in the case of an "im- 
perfect" BGK initial condition where spacetime rapidity 
and the longitudinal flow rapidity are not perfectly cor- 
related. Within the "ideal" Bjorkcn limit the correlation 
is indeed perfect, but its reasonable to expect that devi- 
ations occur. 

The existence of the "ridge" [l9[ provides an experi- 
mental indication for the existence and the size of these 




FIG. 2: (color online) Initial densities in the BGK model (left), as well as the model used in \7) (right). In the BGK case, 
dashed lines represent the rapidity extent of the "excited state" produced by the individual nucleon, while solid lines correspond 
to the cumulative density. See text for model definitions and further explanations 



deviations, independently from the detailed mechanism 
for it's origin: Whatever causes the ridge involves degrees 
of freedom separated in configuration space by ~ fm (the 
jet cone volume) which however are separated in flow ra- 
pidity by the ridge size. (Considerably more than the 
space time rapidity extent of the jet cone). It is then 
likely that such correlations appear in interactions be- 
tween soft (thermalized) degrees of freedom just as they 



appear in the interactions between the system and the 
jet. 

It is reasonable to assume, as an ansatz, the deviations 
are Gaussian, and the density in 77 of matter flowing with 
rapidity y is 



<j>{y, 77) ~ exp 



-(v - vf 



2a% 



(16) 



5 



where a v is a parameter to be determined. Putting in this 
factor instead of the ^-function in Eq. 1131 and integrating 
Eq. \TE\ yields, at 77 = 



1 



2V2vr 



Be 



^V r) (2-2e 2 ^+ (17) 



-vlY 



2tt<t 7 , 



-erf 



Vl 



\/2<7 n 



erf 



VL 



this somewhat unwieldy expression simplifies, at mid ra- 
pidity, to 



(Pz 



cx Be a v /2 o* 



(18) 



The net-vorticity at mid-rapidity then becomes equal, 
roughly to the initial vorticity in the firestreak case multi- 
plied by e rTf >^ 2 a 2 and divided by the longitudinal velocity 
of the nuclcons. 



{p q R ) 


BGK 


(p q R ) 


firestreak 



c(s) 



(19) 



In the limit of — ► the system has no vorticity. While 
only at very low energies (where the formula in Eq. [TS] 
and the BGK picture are untenable as approximations) 
the BGK and Firestreak pictures are comparable, vortic- 
ity at BGK could still be non- negligible provided a n ~ 1. 

It should be underlined that c(s) contains very differ- 
ent physics from a n : In Q, c(s) is interpreted as the 
number of partons into which the energy of the initial 
collision energy is distributed. a v , on the other hand, de- 
pends on the imperfection of "Bjorken" expansion (cor- 
relation between spacetime and flow rapidity). These 
two effects, however, go in the same direction, although 
c(s) ~ (ln-y/s) is much less efficient at diminishing po- 
larization than a small a n . 

Combining c(s) of Eq. QT] with Eq. [ljj] we obtain the 
ratio between BGK and Firestreak expectation, and it's 
dependence on energy and the parameter a^. The result 
is shown in Fig. [4j assuming a v <C ul- The purpose of 
this figure should be taken as an illustration of the sen- 
sitivity of the polarization measure to the longitudinal 
structure of the initial condition, rather than as a pre- 
diction of the polarization in the two models ( as shown 
in [lo| . the small angle approximation used in [?J is in 
any case likely to be inappropriate). As can be seen, the 
effects of c(s) in the Firestreak picture are comparable to 
the effects of a non-negligible a n in the BGK picture only 
at low energies (where the Firestreak picture is thought 
to work better). At top RHIC energy, even at a v of one 
unit, the BGK polarization should be suppressed with re- 
spect to the firestreak expectation with about two orders 
of magnitude. This grows to several orders of magnitude 
for LHC energies. 



Thus, the measurement of the A polarization in the 
reaction plane could be an valuable tool of the initial 
longitudinal geometry of the system. At the moment the 
longitudinal geometry, and in particular how the longitu- 
dinal scale varies with energy (at what energy, and if, and 
how do initial conditions go from "firestreak" to "BGK") 
, is not well understood [20(. This understanding is cru- 
cial for both the determination of the equation of state 
and the viscosity, since longitudinal geometry is corre- 
lated with the initial energy density, and hence to the 
total lifetime of the system and the time in which flow 
observables can form 



21] 



The measurement of the energy and system size de- 
pendence of A polarization in the reaction plane at mid- 
rapidity could be a significant step in qualitatively assess- 
ing the perfection of the fluid, and determining at what 
energy does the system enter fluid-like behavior. 

Connecting the experimental measurement of the A po- 
larization to the initial condition is, however, non-trivial, 
as this observable is sensitive not just to the initial stage 
but also to the subsequent evolution of the system, up to 
the final freeze-out. 

In the next two sections we will qualitatively discuss 
the effect the later stages will have on the final observ- 
able. We will argue that, while the observable is likely to 
be modified by the subsequent evolution, a comparison of 
several kinds of polarization could be useful in obtaining 
information about not only initial conditions, but also 
the mean free path and the freeze-out scenario. 



HYDRODYNAMIC EVOLUTION, 
POLARIZATION AND JETS 



In relativistic hydrodynamics, 



vorticity works some- 



what differently than in the non-relativistic limit [22 
While in non-relativistic ideal hydrodynamics, the con- 
served circulation is defined simply as V x v, relativisti- 
cally the conserved vorticity is 



Q = V x wyv , 



(20) 



where w is the enthalpy per particle. 

In the non-relativistic limit, where w ~ m and 7=1 
the usual limit is recovered. In a relativistic fluid with 
strong pressure and energy density gradients, on the 
other hand, vortices can be created and destroyed even 
in perfectly smooth initial conditions, such as the BGK 
case described in the previous section. 

As vorticity development is a highly non-linear phe- 
nomenon, quantitative details require numerical simula- 
tions. For demonstrative purposes, we include in this 
paper the vorticity which develops when a momentum 
source moving at the speed of light traverses a uniform 
relativistic fluid. This could be an appropriate descrip- 
tion of the thermalized jet energy loss, if the jet loses 
energy fast and locally. The calculation was done using 



6 



BGK 




FIG. 3: (color online) Initial shear in the BK 

a (3+ 1)D hydrodynamical code [HI- The flow vector in 
the x — y — z coordinate system where the fluid is at rest 
(co-moving with the collective flow) is shown in Fig. [5j 
in parallel (left panel) and perpendicular (right panel) to 
the jet direction. 

The simulation shown in Fig. \E\ is based on a jet en- 
ergy loss model that assumes a high momentum gradi- 
ent. This means either considering the jet as a "ball 
of fluid" moving at ultra-relativistic speeds with respect 
to the medium, or as a source injecting a considerable 
amount of momentum into the medium. The second case 
is similar to the source found in [25|, [2(| to lead to the 
formation of Mach Cones. Recent analytical solutions 
within a strongly coupled N=4 supcrsymmctric plasma 
[13] show that these vortex-like structures persist in non- 
equilibrium strongly coupled quantum field theories, sug- 
gesting that, if the QGP at RHIC is strongly coupled, 
they might be found even if the matter surrounding the 



Firestreak 




<p z x x > (b) 



> 

model (a), as well as the model used in [3] (b) 

jet is not quite in local thermal equilibrium. 

It is not surprising that a large initial momentum gra- 
dient, such as that produced by a jet quickly losing en- 
ergy, can introduce vorticity into the system. As shown in 
the included simulation, these vortices are stable enough 
to last throughout the lifetime of the fluid. Perhaps, 
therefore, an interesting polarization measurement to at- 
tempt is to trigger on events with jets, and measure A 
polarization Pj( in the plane perpendicular to the jet pro- 
duction plane. Since vorticity in such events exists inde- 
pendently of the global initial conditions, this measure- 
ment is sensitive only to the mean free path and perhaps 
final state effects. 

Fig. [5] also illustrates how such a measurement could 
be performed: the polarization axis is defined based 
on the jet (high px trigger) direction. Since vortices 
above and below the jet move in opposite directions, in 
a static medium detecting vorticity via polarization mea- 
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FIG. 4: (color online) Ratio of BGK to firestreak predictions as a function of yfs and <j v , the correlation length between 
spacetime and flow rapidity, calculated using Eqn. 1191 and fTTI 



surements would be impossible. 

If, however, the smoke-ring is in a medium undergoing 
transverse or longitudinal expansion, the flow introduces 
a correlation between the A position within the smoke- 
ring (and hence its polarization) and its average momen- 
tum (p) A . Measuring the polarization of moderately high 
momentum but thermal As (~ 700Af eV) in the plane de- 
fined by the A momentum and the jet direction should 
therefore yield a non-zero result. 

We therefore suggest to measure the polarization Pf 
of jet-associated moderate momentum particles, in the 
plane defined by the jet direction and the direction of the 
particle. The observation of this polarization would be 
a strong indication of collective behavior, since it would 
signify jet-induced vorticity. 

Unlike production plane vorticity, jet vorticity does not 
depend on initial condition, but should hold for a wide 
variety of jet energy loss scenarios, provided the coupling 
between the system and the jet, and within the system's 
degrees of freedom, is strong. 

It is not at all clear, however, wether in the strongly 
coupled regime (rather than the perturbative one, on 



which the calculations of [7|] are based) vorticity will read- 
ily transform into quark polarization. The next section 
is devoted to this topic. 

MEAN FREE PATH AND POLARIZATION 

In a perfect fluid angular momentum should go not 
into a locally preferred direction, but into vortices where 
each volume element is locally isotropic in the frame co- 
moving with the flow. Such vortices should imprint final 
observables via longitudinal collective flow ( e.g. odd v n 
coefficients away from mid-rapidity), but not via polar- 
ization. In this regime, the equations derived in the first 
section are no longer tenable since they assume unpo- 
larized incoming particles and a coupling constant small 
enough for perturbative expansion. 

Keeping the first of these assumptions would violate 
detailed balance, while the second assumption is proba- 
bly incompatible with strong collective behavior. Thus, 
if the equilibration between gain and loss terms happens 
instantaneously ("a perfect fluid"), any created polar- 
ization would be instantly destroyed by subsequent re- 



8 



interactions (In the terminology of eq. [H c(s) — > oo, 
through, since the particles would be "infinitely strongly" 
correlated, this would not imply infinite entropy). The 
local isotropy of a perfectly thermalized system was used 
Q to suggest the disappearance of the production plane 
polarization observed in elementary collisions could be a 
signature of deconfmcmcnt. 

A first order correction comes where the size of the ra- 
dius of curvature within the vortex becomes comparable 
to the mean free path l m fp- The anisotropy would then 
be given by the deformation of a volume element of this 
size. Dimensional analysis, together with the insights 
provided in the first section, yields 



(Pg> ~ tanh 



T 



d(pk 
dxj 



(21) 
(22) 



where (pj) is local direction of momentum in the lab 
frame, T is the temperature, and i is any direction where 
a non-zero polarization is expected (production, reaction 
or jet). Note that Eq. |2TI is simply (p X x), taken over 
a homogeneous volume element the size of the mean free 
path. Thus, potentially, the amount of residual polariza- 
tion that survives hydrodynamic evolution (wether from 
initial geometry or from deformation of the system due 
to jets) is directly connected to the system's mean free 
path. 

Measuring the polarization rapidity dependence (in any 
plane, production, reaction, or jet, where it could be ex- 
pected to be produced) could perhaps ascertain the ra- 
pidity domain of the QGP. If the (s)QGP is formed at 
central rapidity, while the peripheral regions are via a 
hadron gas, one should observe a sharp rise in produc- 
tion, reaction and jet plane polarization in the peripheral 
regions 

The problematic aspect of using polarization for such a 
measurement is that it is sensitive to late-stage evolution, 
including hadronization and the interacting hadron gas 
phase. 

As shown in [28|, an unpolarized QGP medium at 
freeze-out will, through hadronic interactions at last scat- 
tering, produce a net production plane polarization due 
to the Hadronic interactions, in a similar way to how un- 
polarized p-p and p-A collisions result in the net hyperon 
polarization. While local detailed balance will inevitably 
cancel out such local polarization, the rather large mean 
free path of an interacting hadron gas, and the consider- 
able pre-existing flow ensure that any interacting hadron 
gas phase should be well away from detailed balance, and 
hence likely to exhibit residual polarization. 

It then follows that the absence of production plane po- 
larization would be a strong indication not only of sQGP 
formation, but of a "sudden" freeze-out where particles 
are emitted directly from the QGP phase. 



The evidence of quark coalescence even at low momen- 
tum [29}, together with sudden freeze-out fits [3(J 31 1, 
make this scenario interesting enough to be investigated 
further using the polarization observable in any plane 
where the vorticity in the hot phase is expected to be 
non-zero (reaction, production and jet). 

If polarization in all directions is consistently mea- 
sured to be zero, including events with jets and within 
high rapidity bins, it would provide strong evidence that 
the mean free path of the system is negligible, and final 
state hadronic interactions are not important enough to 
impact flow observables. A measurement of production 
plane but not reaction plane polarization would provide 
evidence that the initial state of the system is BGK-like, 
and that the interacting hadron gas phase leaves a signif- 
icant imprint on soft observables. (The BGK nature of 
the initial condition can then be further tested by scan- 
ning reaction plane polarization in rapidity). 

The observation of polarization in the jet plane could 
provide a further estimate of the mean free path, and 
show that the jet degrees of freedom are thermalized and 
part of the collective medium. A sudden jump in any 
of these polarizations at a critical rapidity could signal 
a sharp increase in the mean free path, consistent with 
the picture of a mid-rapidity QGP and a longitudinal 
hadronic fragmentation region. Analogously, a drop of 
polarization while scanning in energy and system size 
could signal the critical parameters required for a transi- 
tion from a very viscous hadron gas to a strongly inter- 
acting quark-gluon liquid. 

In conclusion, we have made a few considerations re- 
garding the physics relevant for hyperon polarization 
measurements in heavy ion collisions. We have argued 
that polarization physics is directly connected to sev- 
eral of the more contentious and not understood aspects 
of the system produced in heavy ion collisions, such as 
initial longitudinal geometry and microscopic transport 
properties. We have, however, shown that the polar- 
ization observable can be significantly altered by all of 
the stages of the system's evolution, thereby rendering 
a quantitative description of it problematic. We have 
argued that measuring polarization in several directions 
(reaction plane, production plane and jet plane), and it's 
excitation function and rapidity domain, could shed some 
light on the described ambiguities. 
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FIG. 5: (color online) Vorticity generated by a fast "jet" traversing the system in the positive x direction. The arrows in 
the left panel show the momentum density of fluid elements in the x-y plane, while the contour in the right panel shows the 
x-component of the velocity in the y-z plane. The jet has been traveling for t = 11.52 fm/c through a static medium 26]. The 
dashed arrows in the right panel indicate the expected direction of polarization of the A (out of plane for left panel, tangentially 
in right panel). If the medium undergoing transverse and longitudinal expansion, the A position within the smoke- ring is 
correlated with it's mean momentum. Thus, measuring A polarization in the plane defined by its momentum and the jet 
momentum should yield a positive net result 
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